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Supercooled liquids and dense colloids exhibit anomalous behaviour known as "spatially het- 
erogeneous dynamics" (SHD), which becomes increasingly pronounced with approach to the glass 
transition jTHS] . Recently, SHD has been observed in confined granular packings under slow shear 
near the onset of jamming, bolstering speculation that the two transitions are related [1HH]- Here, we 
report measurements of SHD in a system of air-driven granular beads, as a function of both density 
and effective temperature. On approach to jamming, the dynamics become progressively slower and 
more spatially heterogeneous. The rapid growth of dynamical time and length scales characterizing 
the heterogeneities can be described both by mode-coupling theory [?] and the Vogel-Tammann- 
Fulcher (VTF) equation [8], in analogy with glass-forming liquids. The value of the control variable 
at the VTF transition coincides with point-J [9] [10] the random close-packed jamming density at 
which all motion ceases, indicating analogy with a zero temperature ideal glass transition. Our find- 
ings demonstrate further universality of the jamming concept and provide a significant step forward 
in the quest for a unified theory of "jamming" in disparate systems. 



PACS numbers: 

At low temperature, high density, and low driving, the 
constituent particles in supercooled liquids [1] , dense col- 
loids [3] , and granular packings [IHS] , respectively, are 
nearly locked into a single disordered configuration in 
which motion is heterogeneous in space and time. Dy- 
namics in these systems may be governed by proximity 
to a generic "jamming transition [IT]", beyond which re- 
arrangements cease and the viscosity diverges. Key fea- 
tures of SHD on approach to the transition include un- 
usual correlations |12) in which particles move in one- 
dimensional paths ("strings [13]") that aggregate into 
clusters [H] . and dynamical correlations as measured by 
a dynamic four-point susceptibility \4 [T51H7j . Clusters 
of strings arise naturally in dynamic facilitation [181 US] 
theory and the random first-order theory of glasses [5D] ; 
their shape reflects the fractal nature of dynamical mo- 
tion in these systems [3T]. Strings are also a crucial in- 
gredient in a recent theory of liquid dynamics near the 
glass transition [22] . 

Recent studies demonstrate that close-packed granu- 
lar systems under slow shear exhibit SHD as well fJj-JB], 
bolstering speculation that liquids and granular matter 
share dynamical similarities on approach to the jammed 
state. However, the universality of the jamming hy- 
pothesis has not yet been tested in terms of variation 
in the hallmark dynamical heterogeneities as a function 
of the control parameter. Here, we present the first si- 
multaneous measurements in any experimental system 
of the growth of the cluster correlation length, string 
length, four-point correlation length, and their charac- 
teristic timescales by varying the control parameter. We 
show that the SHD observed in a far-from-equilibrium, 



athermal system of air-fluidized granular beads is essen- 
tially indistinguishable from that observed in thermal 
systems like supercooled liquids and dense colloidal sus- 
pensions. Moreover, we show that theoretical models de- 
veloped for the glass transition can be used to describe 
our granular system, and predict a mode-coupling like 
transition and, more importantly, the jamming transition 
packing fraction, known as point-J [9] [10] from quantities 
characterizing SHD. 

We characterize the spatiotemporally heterogeneous 
nature of dynamics in an athermal, far-from-equilibrium 
system of air-driven steel spheres on approach to jam- 
ming. Compared to sheared or shaken granular systems, 
in which energy is injected at the boundaries, in air- 
driven systems the energy input is uniform in space and 
time. Our granular system consists of a 1:1 bidisperse 
mixture of steel beads of diameters d s = 0.318 cm and 
di = 0.397 cm, with respective masses of 0.130 gm and 
0.266 gm, confined to a circular region of diameter 17.7 
cm. The packing density is varied from an area fraction 
of <f> = 0.597 to 4> = 0.773 by changing the total number 
of beads from 1470 to 1904. Bead motion is restricted 
to rolling within a horizontal plane, and is excited by 
an upward flow of air at a fixed superficial flow speed 
of 545 cm/s. Bead positions are identified by reflecting 
light from their chrome surface to a camera three feet 
above. The duration of experimental runs is 20 minutes. 
By contrast with the molecules in a supercooled liquid, 
here the particles are macroscopic objects driven at ran- 
dom by a continuous input of energy. Consequently the 
speed distributions are non-Maxwellian, and the average 
kinetic energies of the two bead species are unequal. Nev- 
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FIG. 1: (a) Mean square displacement, (b) Number average mobile bead cluster size S c (t). c) Number average string length 
L s (t). (d) Self contribution to the four point susceptibility xts(t)- ( e ) An instantaneous bead configuration where the colouring 
of beads indicates the mobility over a time interval of 12s (the timescale for both maximum cluster size and string length). The 
10% most mobile beads are red; note that they form clusters. Beads moving in strings have vectors superimposed to indicate 
their directional motion. Note that the dynamics are spatially heterogeneous. 



ertheless, as reported previously |23| . the system mimics 
a simple liquid for low 4> and exhibits tell-tale changes in 
the average structure and dynamics at increasing packing 
densities. 

Dynamical characteristics for an example case, <p = 
0.773, are displayed in Fig. [iji-d. The mean-squared 
displacement versus time interval (delay time), Fig. [T^,, 
averaged over all start times and all beads, is perhaps 
the most familiar quantity. It shows ballistic motion at 
short times, diffusive motion at late times, and a plateau 
of sub-diffusive motion at intermediate times that is in- 
dicative of caging. While informative, the mean-squared 
displacement cannot distinguish uniform from heteroge- 
neous dynamics. For this, we perform three measure- 
ments developed for supercooled liquids. The first in- 
volves clusters of beads that are "mobile," i.e. which have 
displacements ranking among the top 10% of all bead dis- 
placements in a given delay time [Mj (in this system, 10% 
gives the largest distinction between mobile and immo- 
bile beads at all packing densities). The configuration 



displayed in Fig. [T£, where beads are colour-coded ac- 
cording to mobility, demonstrates that the mobile beads 
are not distributed at random; rather, mobile and im- 
mobile beads are clustered and spatially separated, indi- 
cating spatially heterogeneous dynamics. One measure 
of SHD is thus the average size of mobile clusters, S c (t), 
defined as the average number of neighbouring mobile 
beads for a given time interval t. The motion within 
a mobile cluster, shown by the displacement vectors in 
Fig. [lj, tends to be correlated into quasi-ld paths called 
"strings [T3] •" Thus a second measure of SHD is the av- 
erage string length, L s (t), defined in terms of the num- 
ber of beads that, within a time interval t, replace the 
initial positions of neighbouring beads to within a toler- 
ance of 0.3d s [15] . Yet a third measure of SHD may be 
constructed from a four-point susceptibility Xi{^)i which 
measures the extent to which the dynamics at any two 
points in space are correlated within a time interval [15j . 
The self-contribution xf S (t) dominates the general re- 
sult [16] and is computed from the variance of the self 
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FIG. 2: (a), (b) Distribution of (a) cluster sizes and 
(b) string lengths for three values of control vari- 
able <j>. The cluster size distribution approaches a 
power law with increasing (j> (solid lines indicate fit 
by power law multiplied by an exponential cutoff), 
while the string length distribution is exponential. 
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overlap order parameter </s(£), which decays from one to 
zero: 

xf(t)=N[( qs {t) 2 )-( qs (tf)}. (1) 

Here, N is the total number of beads, and qs(t) is defined 
as: 

qs(t) = ^w(\r i (t)-r i (0)\)- (2) 

The overlap parameter w is defined by: 

Jl |r,-(t) - r 4 (0)| <Q.U a 
w = < (6) 
[0 |r,-(t) - ri (0)| >0.5d s . 

Here, rj(t) is the position of bead i at time t; averages 
are taken over all beads and over all start times. 

The example results in Fig. [IJ)-d for the cluster size 
S c (t), the string length L s (t), and the four-point suscep- 
tibility Xs«W) a H exhibit well-defined peaks as a function 
of time interval, as found in glass-forming liquids. The lo- 
cations of the peaks indicate the time interval over which 
the dynamics are most heterogeneous, and the heights of 
the peaks indicate the spatial extent or "strength" of the 
heterogeneities. As with glass-forming liquids [HI 1251 - 
|2"T] and colloids [3J, the cluster size and string length are 
largest at the crossover between caged and diffusive mo- 
tion, while XssM ( an d X 4 W) peaks later, in the so-called 
alpha or structural relaxation regime [16j [28] • The ather- 
mal air-fluidized beads therefore exhibit spatially het- 
erogeneous dynamics that is identical to thermal glass- 
forming systems with respect to these three measures. 

Now that spatially heterogeneous dynamics are estab- 
lished for gas-fluidized beads, we turn to their variation 
as a function of control parameter. The distribution of 
cluster sizes at the peak time interval, shown in Fig. [2^, for 
three different packing densities <j>, approaches a power- 
law as <f> is increased. This is consistent with the perco- 
lation of mobile bead clusters; similar power-laws have 



been observed in colloids [3J and in simulations of super- 
cooled liquids [TH [55J near the mode-coupling tem- 
perature. Furthermore, the distribution of string lengths 
at the peak time interval, shown in Fig. [2]d, is exponen- 
tial, , at all values of cj> , where Iq is set by the average 
string length. This is consistent with behaviour reported 
in simulations of several supercooled liquids [131 [551 • 
We note the average cluster size is not much larger than 
the average string length, although the largest clusters 
observed ( 100 particles) are substantially larger than 
the largest string observed ( 30 particles) (not shown). 

Results for S c (i), L s (t), and Xs s W vs ^ are displayed 
in Fig[3ji-c for a sequence of different packing densities 
<j>. When beads are added to the system, the average ef- 
fective temperature also decreases, resulting in a trajec- 
tory in the (<fi ,T e ff) phase diagram that heads towards 
point-J, the zero-temperature jamming transition previ- 
ously found for this system at <fr = 0.83, which is coin- 
cident with the packing density at which the system is 
random close-packed. As the motion becomes more re- 
stricted, the peaks in all three measures of SHD grow 
and move to later times. Therefore, the dynamics not 
only slow down but also become more heterogeneous on 
approach to point-J. Since the SHD functions have ap- 
proximately the same shape when viewed on a log-log 
plot (see data collapse in insets of Fig. |3^,b), this be- 
haviour is fully characterized by the (^-dependence of 
the characteristic or peak time scales {t* Sc ,t* Ls and i* 4 } 

and length scales {CscOsJ, Ci«(*L) and Cx*?*^)}- The 
length (Ls(t*L S ) oc L s is a correlation length for stringlike 
motion, (sc{t*s c ) oc S c is a correlation length of mobile 
particle clusters, and is a correlation length [6J of clus- 
ters of caged particles. On approach to point-J, both the 
characteristic times and the correlation lengths appear 
from Fig. [3] to grow without bound. This is reminis- 
cent of behavior for supercooled liquids as temperature 
is lowered. Though very different, both types of sys- 
tems appear to approach an unusual critical point where 
the growing length scale is purely dynamical, such that 
there is no macroscopic change in instantaneous struc- 
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FIG. 3: (a) Cluster size, (b) string length, and (c) 
self contribution to the four-point susceptibility, all as a 
function of time interval t, for a sequence of area frac- 
tions 0. In order of increasing peak height: = 
0.597, 0.633, 0.647, 0.669, 0.693, 0.722, 0.742, 0.762, 0.773. In- 
sets in (a) and (b) show collapse of all data sets upon scaling 
by the peak heights and peak times. 



Mode-coupling theory, Vogel-Tamann-Fulcher 
4> c =0.79 equation, §o=0.84 
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FIG. 4: (a) Log-log plot of dynamical correlation lengths 
versus 1/|0/0 C — 1| fit with x = 0.40 for £se, x = 0.15 for 
^Ls and x=1.70 for (, X 4. (b) Log-log plot of characteristic 
times t* versus 1/|0/0 C — lj with x — 1.3 for t* Sc , x = 1.1 
for t* Ls , and x = 1.0 for t^ 4 . (c) Semilog plot of dynamical 
correlation lengths versus l/|0/0o — 1| fit with E = 0.10 for 
Csc E = 0.03 for ( Ls , and E = 0.42 for < x4 . (d) Semilog 
plot of all t*s vs 1/1^/^0 - 1| fit with E=0.27. See text for 
fit expressions. Error bars for (a), (c) represent the standard 
error of the measurement. Error bars for (b), (d) represent 
the uncertainty in determining the time at which the function 
is maximum. . 



ture [DIDH32]. 

To further quantify this analogy, the growth of the 
characteristic timescales and dynamical length scales is 
shown in Fig. [4ji-d as a function of packing density. Mo- 
tivated by recent studies [55H55] predicting a power law 
divergence of dynamical lengthscales from mode coupling 
theory [7] (MCT), as well as earlier applications of MCT 
to liquids and colloids, we fit all data to a power-law of 
the form 1/ |0/0 C | X , where both C and x are adjustable 
parameters. As seen in Fig. |4^,,b, excellent fits are ob- 
tained to all data for a single value 0c = 0.79±0.02. This 
value of lies well above the onset of caging and is less 
than the jamming packing fraction, in analogy with well- 
established findings that the mode-coupling temperature 
is below the caging transition but above the glass tran- 
sition temperature [3 , and demonstrates for the first 
time a mode-coupling-like transition in a granular sys- 
tem. In addition to MCT, the glass transition can also 
be well-described by the Vogel-Tammann-Fulcher (VTF) 
equation [S]; therefore, we also fit the characteristic time 
and length scales to the form exp (E/ |0/0o — 1|) , where 



E and O are adjustable parameters. As seen in Fig.|4]:,d, 
excellent fits are obtained to all data for a single value 
0o = 0.84 ± 0.02, consistent with the value of random- 
close packing for the bead system and the value of point- 
J. Since random close packing is the densest random 
packing possible and the point at which all motion ceases, 
a VTF packing fraction of 4>uc P is analogous to an effec- 
tively zero-temperature ideal glass transition, consistent 
with the definition of point-J. This is the first prediction 
of point-J in a granular system from analysis of spatially 
heterogeneous dynamics. 

Our study implies that the behaviour of jammed sys- 
tems, both thermal and athermal alike, may be under- 
stood using the theoretical tools developed for liquids. 
This, in turn, highlights the importance of packing in 
the underlying physics of the glass transition and jam- 
ming. Our results open the door to future theoretical 
insight into the relationship between granular materials 
and supercooled liquids, which might be united by a uni- 
fied theory of jamming. 
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